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We discuss symmetry constraints on the azimuthal oscillations of two-particle correlation (Han- 
bury Brown-Twiss interferometry) radii for non-central collisions between equal spherical nuclei. We 
also propose a new method for correcting in a model-independent way the emission angle dependent 
correlation function for finite event plane resolution and angular binning effects. 
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I. INTRODUCTION 



In the study of heavy ion coUisions two-particle Bose- 
Einstein correlations are an important tool for extracting 
information on the space-time structure of the collision 
zone at freeze-out ||l|. At non-zero impact parameter 
the reaction zone formed by the two overlapping nuclei 
is initially spatially deformed. When viewed along the 
beam direction z, it is longer in the direction y perpen- 
dicular to the reaction plane (defined by the beam axis 
and the impact parameter vector pointing in x direction) 
than in the reaction plane. The reaction plane orienta- 
tion can be determined event-by-event from anisotropics 
in the collective flow of the emitted particles |^, |^ : at 
lower collision energies one exploits the directed flow of 
the protons near projectile and target rapidity ("bounce 
off") while at high energies (where the directed flow 
becomes very weak) one uses the elliptic flow of the pro- 
duced particles at midrapidity 

Having extracted the orientation of the reaction plane 
from the final distribution of the emitted particle mo- 
menta, one can then address the question of their spatial 
distribution relative to the reaction plane by measuring 
10, P two-particle correlations as a function of the az- 
imuthal emission angle $ (i.e. the direction of the trans- 
verse momentum vector K±^ of the emitted particle pairs 
relative to the impact parameter 6) |l^, |ll|. Com- 
plementing the spectral information on the momentum- 
space structure of the source with space-time information 
from the correlation functions severely constrains models 
for the dynamical evolution of the reaction zone . For 
noncentral collisions interesting questions which can be 
addressed in this way are the origin and manifestation 
of anisotropic collective flow and its consequences for the 
space-time evolution of the fireball, from which informa- 
tion about the intensity of rescattering effects and the 
degree of thermalization in particular during the early 
stages of the colHsion can be extracted (see, e.g., [Q). 

For quantum statistical correlations between, say, iden- 
tical pions the measured two-pion correlation function 
C{q,K) is related to the emission function (single-pion 



phase-space distribution at freeze-out) S{x, K) by j|] 



C{q,K) = l+ (e*'J-(^-^*)\ (K 



(1) 



Here q=Pi—p2 and K = ^{pi-\-p2) are the relative 
and average momentum of the pion pair, respectively, 
f^—{Pi+P2)/{Ei-\-E2) is the velocity of the pair, and 
the average (...) is taken with the emission function: 



(/(^)>W- 



/ d'^xfix)S{x,K) 
J d^xS{x,K) 



(2) 



If the space-time structure of >5'(a;, K) can be approx- 
imated by a Gaussian, the resulting correlator C{q,K) 
is again a Gaussian in the relative momentum q and can 
be fully characterized by six HBT radius parameters R^j 
which are functions of the pair momentum K ||l| : 



K)^l + exp 



E 

i,j—o^s,l 



(3) 



Here q is conventionally decomposed into orthogonal 
components along the beam direction [l — longitudinal), 
parallel to the transverse pair momentum K^^ (o = out) 
and along the remaining third direction {s = side). In 
this (osZ)-frame the pair velocity has components (3 — 
{f3±,0,(3i). The radius parameters R'^j are then calcula- 
ble from the spatial correlation tensor 

V(-f^) = {x^x,){K) - {x^,){K){x,)iK) = (i^i,), 

(4) 

(/i, u = 0,1, 2, 3), which describes, for pairs with momen- 
tum K , the widths in space-time of the emission function 
S{x,K) around the point of highest emissivity |l|. The 
spatial correlation tensor is specified in coordinates 
attached to the reaction plane: x^ = z is the beam di- 
rection, xi= X is the direction of the impact parameter 
b, and X2=y points perpendicular to the reaction plane. 
For the spatial correlation tensor this choice of coordi- 
nates is natural since the reaction plane is a symmetry 
plane for the collision. The relations between i??- and 
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are § 

Rl - i (511+^22) (^11 -5'22)cos(2$)-5i2sin(2$) 
^ i(5n+522) + 5(^11-^22) cos(2$) + 5i2sin(2$) 

-2/?_l(5'oicos$+S'o2 sin $) 5*00, 
= 5i2COs(2$)-i(S'n-522)sin(2$) 
+/3_l(S'oi sin $-5*02 cos<I>), 
-R? = Sza, — 2/3; 5o3 + PfSoo, 
RIi = 



(523-A -5*02) sin $ 
" " ~ ^ (Sia-AS-oOsin^. 



^o; = (5i3-A5oi)cos$ 

Rll - (523-A5o2)cOS$- 



(5) 



The radius parameters R^^ are functions of the pair 
rapidity Y = i ln[(l+/3i)/(l— /?;)], the magnitude K^^ of 
the transverse pair momentum, and its angle relative 
to the reaction plane (the azimuthal emission angle). In 
Eqs. (^) we only indicated the explicit $ dependence aris- 
ing from the azimuthal rotation of the {osl) system rela- 
tive to the reaction-plane-fixed (xyz) system: 

Xo = X cos $ -f ?/ sin $, = — xsin<i> -I- y cos$. (6) 

In addition to this explicit $ dependence, there is an im- 
plicit one pO| arising from the dependence of the emission 
function S{x, K) — S{x, y, z, t; Y, K±, $) on the emission 
angle this generates a $ dependence of the compo- 
nents of the spatial correlation tensor 5*^,^. In this note 
we work out the symmetry constraints on the $ depen- 
dence of Sfiu and study their implications on the $ de- 
pendence of the HBT radius parameters after the explicit 
$ dependence shown in Eqs. (||) is folded in. 



II. SYMMETRIES OF THE EMISSION 
FUNCTION AND SPATIAL CORRELATION 
TENSOR 

For spherical colliding nuclei the emission function is 
symmetric under reflection at the reaction plane: 

I : S{x, y, z, t; Y, K^, $) === S{x, -y, z, t; Y, K\, -$). 

(7) 

This leads to the following symmetry relations for the 
spatial correlation tensor: 



with 



(-1) 



(8) 



(9) 



Thus, symmetry I relates the components of 5^1, at emis- 
sion angle $ with those at angle —3' at the same pair 
rapidity Y and transverse momentum K±. S02, S12 and 
523 are odd under this symmetry (6*1 = — 1), all other 
components are even (9i = -t-1). 



If the two nuclei have equal mass, the emission function 
is also symmetric under interchange of projectile and tar- 
get. In the center of mass system centered at the collision 
point, this translates into a point reflection symmetry at 
the origin: 

II: S{x,y,z,t;Y,K^,^) 

^ S{-x,-y,^z,t;-Y,K^,'i> + n). (10) 

For the spatial correlation tensor this implies 

II : S^,{Y, K^, $) =. 92 S^,i-Y, K^,^+7r), (11) 

with 



92 = (-1)''"°^ 



(12) 



Symmetry II relates Sfj,i, at emission angle $ for forward- 
going pairs (y > 0) with at emission angle ^ + tt 
for backward-going pairs (Y < 0), and vice versa. For 
midrapidity pairs {Y = 0) it relates the spatial correlation 
tensor at emission angles $ and $-|-7r, providing a use- 
ful second constraint on the emission angle dependence. 
Sqi, Sq2 and 5*03 are odd under this symmetry {02 = — 1) 
while all other components of 5'^,^ are even (6*2 = +1). 

Finally, at very high collision energies the source 
is expected to be approximately invariant under lon- 
gitudinal boosts within an extended rapidity interval 
around Y = 0. If this is the case, the emission func- 
tion S{x,K), when expressed in terms of longitudi- 
nal proper time r = \^t^—z-^ and space-time rapidity 
ri = ^\ii[{t+z)/{t—z)], depends only on the difference 
rj—Y between the space-time and momentum-space ra- 
pidities. For equal projectile and target nuclei it then 
must be an even function of r]—Y^ i.e. invariant under a 
simultaneous sign change of Y and rj. With z^rsinhr; 
and t = T cosh rj this implies 



III : S{x, y, z, t; F, K^, $) = 5(x, y, -z, t; -Y, K^, $) 



and 



III : S^.{Y, K^, $) = 03 5^.(-r, Xi, $). 



with 



= (-1) 



(13) 



(14) 



(15) 



Combining symmetries II and III allows to relate the 
spatial correlation tensor at angles $ and <I>+7r for all 
rapidities Y . For boost-invariant sources, the terms with 
03 = —1 (i.e. 5*03, 5*13, and 523) vanish at 1" = 0. We note 
that the symmetry ( |l3|) also applies to sources without 
boost-invariance if they exhibit spatial and momentum 
anisotropics (i.e. ^-dependence) already at zero impact 
parameter, such as fully central collisions (no spectators) 
between deformed nuclei (e.g. U-l-U). In this case the 
source is symmetric under the simultaneous reflection 
of coordinates and momenta at the transverse plane at 
z = 0, in agreement with Eq. (13) p^. 
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TABLE I: Consequences of symmetries I and II (see text) 
for the azimuthal Fourier expansion of the spatial correlation 
tensor S^i, at midrapidity Y = Q. The last column lists the an- 
gles $ in the first quadrant where Si_i,u{Y — 0, K±, "I>) vanishes. 
The notation follows the one introduced in pXh. 
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TABLE II: Consequences of symmetries I, II and III (see text) 
for the azimuthal Fourier expansion of the spatial correlation 
tensor S^i, for boost-invariant sources. The last column lists 
the rapidities Y and angles "I> in the first quadrant where 
Si_ii, (Y, K± , 3>) vanishes. 
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We will concentrate here on the consequences of the 
combination of symmetries I and II at Y = and of the 
combination of all three symmetries at any Y. The for- 
mer case is relevant for two-pion correlations at midra- 
pidity in low-energy collisions between equal spherical 
nuclei, the latter case applies to high energy collisions, 
such as those studied at the heavy-ion colliders RHIC 
and LHC. Symmetry I alone is less restrictive and is the 
only useful one when significantly away from midrapidity 
(in particular in the projectile and target fragmentation 
regions). 



III. AZIMUTHAL FOURIER DECOMPOSITION 
OF THE SPATIAL CORRELATION TENSOR 

The above symmetries constrain the $ dependence of 
the components of the spatial correlation tensor (and 
thereby the implicit $ dependence of the HBT radius 
parameters in Eq. (|^)). Correspondingly, certain expan- 
sion coefficients will vanish in an azimuthal Fourier ex- 
pansion of S'p^. Let us write generically S{^) for the 
^-dependence of a given component S'^jy. Being a real 
function it has the following Fourier decomposition: 



S{^) = Co + 2 ^ [C„ cos(n$) 5„ sin(n$)] , (16) 

r - 



n=l 

— S{^) cos(n$) , 



2ti 



S'($) sin(n$) , 



Symmetry I implies 



c 



- for all n, 
for all n. 



(17) 



(18) 



At Y = Q, symmetry II eliminates even or odd terms in 
the Fourier series: 

6*2 — +1 =^ C„, 5„ = for odd n, 

62 = —1 C„, 5„ = for even n. (19) 

In fact, Eq. ([ill ) implies a stronger result: 

9-2 = +1 Cn,Sn are odd (even) functions of Y 

for odd (even) values of n, 

02 — —1 =^ Cn,Sn are odd (even) functions of Y 

for even (odd) values of n. (20) 

This will be used in Section 0. At F = Eq. (|l9|) follows 
from (1^). 

Table Q lists the Fourier expansions for the components 
of the spatial correlation tensor at midrapidity which re- 
sult from the combination of these two symmetries. Fol- 
lowing ||l^ we have used the fact that Sn and S22 have 
structurally identical Fourier expansions and combined 
them into A — ^{811+822) and -B = ^(S'li— S'22), which 
are the combinations entering in Eqs. (a). 

For boost-invariant sources, we can combine symme- 
tries II and III (Eqs.(|0|) and (|l|)) to obtain 



(21) 



V(y, K^, $) = 02^3 S^AY, K^, $ + tt). 
For the Fourier coefficients this implies 



72173 



-I C„, Sn = for odd n, 

-1 =^ Cn.Sn = for even n. (22) 



In contrast to Eqs. (^^ this is now true for all rapidi- 
ties Y. The corresponding Fourier expansions are listed 
in Table ^ Note that according to Eq. (|2^) the non- 
vanishing coefficients Gn, Hn and /„ in Table || are odd 
functions of rapidity and thus vanish at y = 2. 
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IV. FOURIER EXPANSION OF THE RADIUS 
PARAMETERS 



We will now combine the above implicit $ dependence 
of the spatial correlation tensor with the explicit $ de- 
pendence shown in Eqs. (^. When studying the combi- 
nation of symmetries I and II for sources which are not 
invariant under longitudinal boosts, we must restrict our 
attention to midrapidity pairs with Pi —0. For simplicity, 
we also set (3i—0 in the boost-invariant case; this means 
that we are studying the correlation radii in the longitudi- 
nally comoving system (LCMS) General expressions 
for /3; ^ are easily obtained by boosting from the LCMS 
to a fixed longitudinal reference frame and can be found 
in Refs. |l^. For boost-invariant sources i?^; = in- 
dependent of rapidity |l^; without boost invariance 
this is generally not even true at midrapidity (see Eqs. (|^) 
and Table ||) . As shown in Ref . , a non-zero value for 
R^i arises naturally if the longitudinal major axis of the 
source ellipsoid is tilted away from the beam direction; 
longitudinal boost invariance forbids such a tilt. 

Using the symbols introduced in Tables | and || and 
setting Pi =0, Eqs. (||) simplify to 



Rn 



Rqs = 

Rj = 

Roi = 

R^, — 



A-Scos(2$)-Csin(2$), 

A + 5cos(2$) + Csin(2$) 

-2P±{E cos^+F sin^) + PlD, 

Ccos(2$) - Ssin(2$) -K/3_L(£^sin$-Fcos<I>), 

J, 

cos $ + / sin ^ - P±G, 
/cos* - i/sin$. (23) 



Here A, B, . . . , J are functions of $ whose Fourier expan- 
sions are given in Tables || and ||. For a boost-invariant 
source G($), -ff($) and /($) (i.e. .R^, and i?^;) vanish at 
Pi = Y = 0. A comparison of Tables g and || shows that at 
1^ = all other S'^i^ components entering Eqs. ( p3| ) have 
exactly the same Fourier expansion with and without lon- 
gitudinal boost invariance. We may therefore investigate 
Eqs. ( p3|) on the basis of the expansions listed in Table || 
and recover the boost-invariant case later by simply set- 
ting R^i = R^i=0. 

In Ref. pl| we studied the limit of vanishing implicit 
$ dependence (i.e. Sfj_i, does not depend on $). Table | 
shows that in this limit only the diagonal elements and 
Si3 (i.e. A, B, D, H and J) are non-zero. As discussed 
in , this limit requires only that space- momentum cor- 
relations (e.g. due to collective flow effects) are weak. A 
more general situation was analyzed in [|lO| where terms 
up to rt = 2 were kept in the Fourier expansion of the 
HBT radii, but in that work the $ dependences of the 
emission duration and time-space correlations were ne- 
glected relative to those of the spatial components Sij. 
We here remove both of these approximations. 



Inserting the expansions in Table | into (^3|) and using 

cosn$cosm$ = i [cos(?i— to)$ + cos(?i-|-m)<I>] , 
sinn<I> sinm<I> — i [cos(n— to)$ — cos(rt+m)<I>] , 
cosn<I> sinm<I> = ^ [sin(n-|-m)<I> — sin(n— to)$] ,(24) 

we see that at midrapidity the HBT radius parameters 
_R^(y=0, K^, $) have the following Fourier expansions: 



Rl 


- Rlfi - 


^2E„= 


=2,4.6,. 


..^Lcos(n$), 


Rl 






=2,4,6,. 


i?2„cos(n$), 


n2 






=2,4,6,. 


,.^L,nSin(n$), 


Rl 


= RIo ^ 




=2,4,6,.. 


Rln cos(n$), 


Rli 






= 1,3,5,. 


,.^o;,ncos(n$), 


Rli 






= 1,3,5,. 


,.-Rl«sin(n$). 



Due to the symmetries of the emission function, the HBT 
radius parameters are sums of either cosine or sine terms, 
involving either even or odd multiples of the emission 
angle $, but no mixtures of different such terms. As 
a consequence, Rl^ vanishes at both $ = 0° and 90° (i.e. 
its leading contribution features a second order harmonic 
oscillation as a function of the emission angle). Rli and 
RIi in general exhibit leading first order harmonic oscil- 
lations with zeroes at 90° and 0°, respectively. For 
a boost invariant source they vanish identically. 

The Fourier coefficients i?„ „ are functions of K±. We 
now list them up to order n = 2. The <f>-independent 
terms are given by 



s.O 

.2 

■0,0 



Rf.Q 



Aq — B2 — C2, 

Ao + 52 + ^2-2/3^(^1 
Jo- 



PlDo, 



(26) 



The coefficients of the first order harmonics are 



Roi,i — ^{Ho + H2 + I2) — P±Gi, 
RliA = \{-Hq + H2 + h)- 



(27) 



For a boost-invariant source these vanish. The term ^ Hq 
describes the tilt of the emission region relative to the 
beam axis which was discussed in . The second order 
harmonic oscillations have amplitudes 



R 



S.2 



rI ■ 



R' 

R 



OS, 2 
2 

7,2 



= A2-^{Bo + Bi + Ci), 
= A2 + ^{Bo + Bi + Ci) 

-p^{Ei +E3-F1+ F3) + PID2, 

- li~Bo + Bi + Ci) 

= J2. 



/9. 



i (El — E3 — Fi 



F3), 
(28) 



If the emission duration D = {t^) is independent of emis- 
sion angle (152,4.6,... ~0) and all higher order harmonics 
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n > 3 of the spatial correlation tensor are small, these 
amplitudes fulfill the approximate "sum rule" ||l0|| 

p2 p2 1 o p2 

-"■0,2 " ^s,2 + ^-ftos,2 — 

2(^4+^4) - 2^^{E3+F3) + (3lD2. « 0, (29) 

Note that the leading first order harmonics of (tx) and 
(ty) , which describe how the transverse positions are cor- 
related with time at freeze-out, cancel in this "sum rule" . 
If the data satisfy this "sum rule" for all values of K± 
resp. /3_L, one may conclude (barring unlikely acciden- 
tal cancellations among the terms) that D2, E3, F3, B4 
and C4 all vanish. In this case the azimuthal oscillation 
amplitudes of the transverse HBT radii reduce to 

^1,2 = M - ^Bo, 
Rla = A2 + ^Bo-l3^{Ei-Fi), 
Rls,2 = -\B^ + \P±{Ei~F^). (30) 

The term ^{Ei—Fi) is the leading (first harmonic) 
contribution to the correlation {{x—y)t) between emis- 
sion points and times. In a hydrodynamic model this 
term reflects the geometric manifestation of elliptic flow, 
namely that the freeze-out radius increases with time 
more rapidly in x than in y direction ||l6|| . Since it comes 
with an explicit factor of one may be able to isolate 
it using the i^_L -dependence of the azimuthal oscillation 
amplitudes ( pO| ) at small K±^. 

Finally, in the absence of dynamical space-momentum 
correlations, all implicit $ dependences (i.e. all higher 
harmonics in Table |) are expected to vanish, leading to 
the "geometric relations" [|lO| 



R2 


^0, 




r2 

-"-0,0 






Ruo = 


Jo, 






- ~Rsl,l = 




Rl,2 = 


-Rl,2 = " 





In this case all five non- vanishing components of the spa- 
tial correlation tensor can be separated . 



V. CONSIDERATIONS FOR A FINITE 
SYMMETRIC WINDOW AROUND y = 

The results quoted so far were derived at midrapid- 
ity y = since, at least in the absence of longitudinal 
boost invariance, symmetry II can only there be used 
to constrain the azimuthal Fourier expansion of S^^, by 
eliminating either even or odd terms in the sums over 
n (see Eq. (p^)). In practice, statistical hmitations ren- 
der strict cuts on the pair rapidity Y quite painful. It 
is therefore important to assess the necessary modifica- 
tions if the data are collected in a finite size rapidity in- 
terval around 1^ = 0. We now prove the important result 



that, as long as the HBT radii are obtained from av- 
eraging over a symmetric rapidity interval around Y — Q, 
the general form ( p5| ) of their Fourier expansions remains 
unchanged. On the other hand, equations ( p6|) - (^8|) re- 
ceive additional contributions which, at leading order in 
the width Ay of the rapidity interval, grow quadratically 
as (Ay)^; this can be used to eliminate them by varying 
Ay and extrapolating to Ay = 0. 

As noted in Eq. (pO[), the point reflection symmetry 
(|To|), ( p] ) allows to classify the Fourier expansion coeffi- 
cients of the spatial correlation tensor as either even 
or odd functions of rapidity Y . The odd terms vanish at 
y = 0, but do not do so any longer at y^O. However, 
when calculating the HBT radii from S^i, according to 
Eqs. (|^) and averaging them over a finite, but symmet- 
ric rapidity interval around y = 0, terms which are odd 
in y average to zero. Therefore, there are no new con- 
tributions to i?^, i?Q and E?gg in this case. i?^; and 
i?^;, on the other hand, contain at y 7^ additional terms 
beyond those listed in Eqs. ( p3| ) which are multiplied by 
either one or two powers of /?/ . When multiplying an odd 
Fourier coefficient by /9/, the result is even in f3i (respec- 
tively y) and does not average to zero across the rapidity 
interval Ay. In fact, at leading order in Ay, its average 
is ^ {0f) which grows quadratically with Ay. 

Let us now look at how these extra terms modify 
the Fourier expansions given in the last three lines of 
Eq. (H). We begin with i?f = 5*33 - 2/3/5o3 + PfSm and 
average it over the symmetric interval Ay. Table ^ tells 
us that the Fourier coefficients of S'33 with odd values of 
n are odd functions of Y and thus average to zero. For 
S'03 the coefficients with odd n are even functions of y, 
but since S'03 is multiplied by /3/ , these odd n terms again 
average to zero. The same is true for the last term where 
the factor 0f preserves the y-refiection symmetries of the 
expansion coefficients. Altogether, the rapidity-averaged 
longitudinal radius ) continues to have only even n 
terms in its Fourier expansion, just as Eq. ( p5| ) states 
for Y — {). In the same fashion one also shows that the 
rapidity- averaged radius parameters {R^i) and (i??^ con- 
tinue to have the same Fourier expansions as in (Pq). In 
other words, averaging the HBT radii over a finite, sym- 
metric rapidity interval around Y — Q preserves the gen- 
eral structure ( p5| ) of their azimuthal Fourier expansions. 

When expressing the Fourier components of the rapi- 
dity-averaged HBT radii in terms of the harmonic coef- 
ficients of S"^!,, new terms arise, and Eqs. (p6|)-(p8|) are 
modified. We only list those equations whose structure 
changes: 

(iJ^o) = (Jo)-2(AGo) + (A'^o), 

{RI2) - (J2)-2(AG2) + (A'i?2), 
{Rli,i) - \{Ho+H2+h~ Pi{Eo+E2+F2)) 

-l3^{Gi-l3iD{), 
{Rls) - \{-Ha+H2+l2~ l3i{-Eo+E2+F2)). (32) 

Here the angular brackets denote the average over the 
symmetric rapidity interval Ay. All terms involving one 
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or two explicit factors (3i vanish quadratically as AK - 
in which limit Eqs. ( pq) - ( |2^ ) are recovered. 



VI. WHAT IF THE SIGN OF THE IMPACT 
PARAMETER CAN NOT BE DETERMINED? 

If the orientation of the reaction plane is reconstructed 
from an even Fourier component of the single particle dis- 
tribution (e.g. from the elliptic flow coefficient V2 as is 
the case at RHIC), the direction of the impact param- 
eter vectoir b has a sign ambiguity, i.e. after aligning 
events according to their reaction plane the event sample 
contains equal contributions from collisions with impact 
parameters b and —6. This ambiguity does not exist if 
the event plane is reconstructed from the directed flow 
coefficient vi (as one does at the AGS and SPS) whose 
sign has a one-to-one correlation with the direction of b 
within the reaction plane. 

If events with impact parameters b and —b are equally 
mixed, the effective source function is symmetric under 
the exchange b —b which is equivalent to an azimuthal 
rotation by 180°: 

Ila : S{x,y,z,t;Y,Kj^,^) 

= S{-x,^y,z,t;Y,K^,'i> + Tr). (33) 



For the spatial correlation tensor this implies 



Ila : Sf^AY, K^,<i>)= O^a S^,,{Y, X^, $+7r) 



with 



(-1) 



(34) 



(35) 



One easily checks that the sign 6*20 is, in fact, equal to the 
product 02^3 of the signs under symmetries I and III, as 
tabulated in Table ||. Correspondingly, the general form 
of the Fourier expansions of 5*^,^ and of the HBT radii 
are exactly the same as those listed in and resulting from 
Table || for a longitudinally boost-invariant source. We 
see in particular that a sign ambiguity for the direction 
of the impact parameter automatically leads to vanishing 
cross terms R^i and i?^;. 



Since the reaction plane is reconstructed event-by- 
event from the anisotropics of the single particle momen- 
tum distribution |2[ ||, ^, its orientation is only known 
with a finite statistical accuracy controlled by the num- 
ber of particles used in the reconstruction process. Corre- 
spondingly, in a statistical average over the event sample 
the true reaction plane angle 7/;^' is distributed around 
the reconstructed one ipm by a probability distribution 
(see Eq. (9) in [|) 



Pii'm-lpR) 



dP 



dv' 



X exp 



d{m{->pm-ipR)) J 27r(T2 

2VmVin cos (m(V'm-'0i?)) 



2a' 



,(36) 



with width cr^ = {w^)/{2M{w)'') (where M is the num- 
ber of particles per event and w is an arbitrary weight 
function (e.g. w = 1 or w=p±) used in the analysis), m 
denotes the order of the Fourier component of the single- 
particle spectrum used to extract the reaction plane, and 
Vm is the corresponding Fourier coefficient; the cases 
TO = 1 (directed flow) and m = 2 (elliptic flow) are rel- 
evant in practice. Correspondingly, a measurement of 
the two-particle distributions N{q,K) and D{q,K) at 
fixed emission angle <& = ^—^pm relative to the recon- 
structed event plane corresponds to an average of the 
real two-particle distributions over a range of emission 
angles ^—^r relative to the true reaction plane, where 
the average is taken with the distribution (^6|). The av- 
eraging reduces the azimuthal dependence of the corre- 
lation function (and of the HBT radii extracted from it) 
and must be corrected for, before comparing to models. 

An additional smearing which goes in the same direc- 
tion arises from the binning of the data in By sum- 
ming the data over all emission angles (f) within an in- 
terval of width A centered at <&, one effectively performs 
an additional smearing of N{q, K) and D{q, K) over the 
azimuthal emission angle with the distribution 



/A(0-<i>) 



1 

A^ 



6-$ + iA) (iA-0+$) . (37) 



VII. CORRECTIONS FOR BINNING AND 
FINITE EVENT PLANE RESOLUTION 



The two effects can be combined by folding the distri- 
butions p and /a and averaging the true correlated and 
mixed pair distributions with 



Experimentally, the two-pion correlation function is 
obtained as the ratio of correlated pairs, N{q,K), 
and uncorrelated (mixed event) pairs, D{q,K). In an 
azimuthally-sensitive analysis, one constructs these dis- 
tributions for a given selection on emission angle How- 
ever, finite binning in $ and uncertainty in the experi- 
mental estimation of the reaction plane tend to dampen 
the azimuthal dependencies of the observed ( "raw" ) dis- 
tributions Ncxp{q,K) and Doxp(g, -K^)- In this Section 
we present a model-independent procedure to correct for 
these effects. 



dip pW /a(0- V') 



1 

A 



A/2 

d9 p{4)~9). 

-A/2 

(38) 



The above azimuthal averaging affects the numerator 
N and denominator D separately. Suppressing the de- 
pendence on K± and Y for clarity, the measured angular 
dependence of the correlated pairs relative to the recon- 
structed reaction plane ipm, -^oxp(Q': ^—'4'm)i is related to 
their true angular dependence relative to the real reac- 
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tion plane i/jr, iV(qf, ^-V'fi), by 

iVcxp(g, $--0™) = f d<p N^iqA-'^R) 

J —TT 



where 



X p(($-V„O-(0-^^)), (39) 



„0_^^+A/2 

NA{q,cl)-^R)^ / N{q,B)de (40) 

"'0-iAh-A/2 



J 



denotes the effect of summing the data in angular bins 
of width A. An analogous pair of equations holds for the 
measured and true uncorrelated pairs in the denomina- 
tor, £)cxp(9, ^-V'm) and D{q,(j)-il;R). 

The task at hand is to extract the true angular de- 
pendence on ^—ipn from the measured dependence on 
^j-'i'm where j labels the angular bins centered at an- 
gles $j relative to the reconstructed reaction plane. To 
this end we Fourier decompose the measured quantities 
N and D for each value of q. For example 



7Vexp(g, $-^™) = A^o^P(9) + 2 J2 KTl?) cos(n(<i>-V'„)) + N^^fiq) sin(n(<i>-^„))] , 

-, "bin 

N'^'^Piq) = (7V,,p(g, $) cos(n<i>)) = ^ iVexp(g, $,) cos(n<i>,), 

^biii . T 
J = l 

^ "bin 

7V|,^„P(q) = (iVexp(g, $) sin(n$)) = — ^ N.^^iq, $,) sin(n$,), (41) 



nhi, 



where ribin denotes the number of angular bins (for finite ribin only Fourier components with n < nbin are meaning- 
ful). We further imagine doing the same for the corresponding true and binned quantities corrected for event-plane 
resolution: 

"bin 

iV(q, ^-iJjr) = Noiq) + 2 [NcAq) cos(n($-V'fl)) + N.^q) sin(n($-V'i?))] , 



N^iq, 0) = <(q) + 2 [<„(q) cosinO) + N^^q) sin(n0)] . 

I 



(42) 



n=l 



Analogous expressions hold for the mixed pairs in the 
denominator D. Inserting the Fourier expansions ( p^ ) 
into Eqs. ( |39| ) and (^0|), using that the distributions p and 
/a are even functions of their arguments, and comparing 
the result with Eq. (^ij) one easily finds for all n and both 
series of coefficients (a = c or s) 



KJiq) 



Na.niq) 



sin(nA/2) 



nA/2 ' 
Knil) (cos {n{ijm-iJn)))p- 



(43) 



The factors (cos {n(iprn—4'R)))p: arising from an average 
over the event plane distribution ( |3^ ) , are the well-known 
correction factors for event-plane resolution arising in the 
process of extracting the anisotropic flow coefficients u„ 
from the single particle spectrum |^ ||, ^. 

With the results (|43|), the numerator N and denomina- 
tor D of the correlation function at each measured angle 
<i>j and relative momentum q can now be easily corrected 
for the effects of angular binning and finite event plane 
resolution (setting again ipm = 0): 



iV(g, $,) = A^exp(q, $,) + 2 ^ Cn,™(A) [N^'^fiq) cos(n<i>,) + N^^^iq) sin(n<i>,)] , (44) 

I 

with correction parameters Cn,m(A) given by the simple expression 

n^/2 , , 

sm(nA/2)(cos(n(V'm-V'ii)))p 
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A similar equation holds for the uncorrelated pairs in the 
denominator D. Since the right hand side of Eq. ( ^ ) 
involves only experimentally known quantities, the cor- 
rection algorithm is model independent. The sums over 
n go over all allowed values; if m is even (i.e. the sign of 
the impact parameter is not known), both N and D are 
symmetric under azimuthal rotations by 180° and only 
even values of n are summed over. Contrary to the HBT 
radius parameters or to single-particle flow measures [Q , 
N and D have no unique symmetry under <!>—!■—$; thus 
in general both sine and cosine terms contribute to the 
sum in Eq. (||). 

After applying the algorithm ( |44| ) to the data, the ra- 
tio C(q, K) = N{q, K)/D{q, K) gives the corrected two- 
particle correlation function from which all angular bin- 
ning and event plane resultion effects have been removed. 
The true emission angle dependence of the HBT radius 
parameters can thus be directly extracted from a Gaus- 
sian fit with Eq. (||) to this function C(<7, K). 



sources they also apply at Y ^Q. The structure of these 
expansions is preserved if the data are averaged over a 
symmetric finite rapidity interval around y = 0. They 
provide a basis for fitting the azimuthal emission an- 
gle dependence of experimentally determined correlation 
functions. A model- independent correction of the mea- 
sured two-pion correlation function for event plane reso- 
lution and angular binning effects is possible and given 
in Section^ Equations @"(||), (|3|) and @ relate 



the oscillation amplitudes extracted from the thus cor- 
rected correlation function to the leading harmonic coef- 
ficients of the spatial correlation tensor and allow to con- 
strain models for the emission function using azimuthally 
sensitive HBT data. Under favorable conditions spatial 
and temporal aspects of the emission function can be 
separated. 
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